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ABSTRACT 


This report is on the third phase of the development of the computer codes for 
scattering by coated bodies that has been part of an ongoing effort in the 
Electromagnetics Laboratory of the Electrical Engineering and Computer Science 
Department at the University of Illinois at Chicago. The work reported here discusses the 
analytical and numerical results for the scattering of an obliquely incident plane wave by 
impedance bodies of revolution with phi variation of the surface impedance. 

Integral equation formulation of the problem is considered. All three types of 
integral equations, electric field, magnetic field and combined field, are considered. 
These equations are solved numerically via the method of moments with parametric 
elements. Both TE and TM polarization of the incident plane wave are considered. The 
surface impedance is allowed to vary along both the profile of the scatterer and in the phi 
direction. 

Computer code developed for this purpose determines the electric surface 
current as well as the bistatic radar cross section. The results obtained with this code have 
been validated by comparing the results with available results for specific scatterers such 
as the perfectly conducting sphere. Results for the cone-sphere and cone-cylinder-sphere 
for the case of an axially incident plane have been validated by comparing the results 
with the results with those obtained in the first phase of this project. In this report results 
for body of revolution scatterers with an abrupt change in the surface impedance along 
the both the profile of the scatterer and the phi direction are presented. 
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1. INTRODUCTION 


The availability of new artificial dielectric and absorbing materials during the 
last decade has resulted in increased activity and interest in the formulation and analysis 
of reflection, transmission and scattering by such materials and bodies coated with these 
anisotropic materials [Graglia et. al., 1984, 1987a, 1987b] The very nature of these 
problems require that, except for very special cases of material parameters and 
geometries, a numerical solution be performed. Finite difference, finite element and 
integral equation formulations are the techniques employed most frequently in such cases. 
In this study the scattering problem is formulated in terms of an integral equation. This 
technique well suited to problems involving anisotropic and inhomogeneous materials. 
The integral equation formulation leads to a set of coupled integro-differential equations 
which can be solved numerically using the method of moments [Harrington, 1968]. From 
a computation stand point a method of moments solution of an integral equation is 
particularly effective when parametric elements are employed [Graglia, 1988a, 1989]. 

This report describes the details of the development of a set of computer codes 
for scattering by a body of revolution with arbitrary surface impedance. This is the third 
phase of an ongoing effort at the University of Illinois at Chicago to develop codes that 
would allow the numerical computation of electromagnetic scattering by and the radar 
cross section (RCS) of an arbitrary body coated with a variety of artificial materials. The 
present phase of the study presented in this report is the development of analytical results, 
numerical codes and results for the scattering of an obliquely incident plane wave by a 
body of revolution with an impedance boundary condition. This phase differs from the 
previous phase in that the surface impedance is allowed to vary not only along the profile 
of the scatierer but also in the phi direction. 
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The analytical results presented in Chapter 2 are applicable to a scatterer which 
has a body of revolution symmetry. The profile of the scatterer is readily described by a 
generating curve described parametrically in cylindrical coordinates by a pair of functions 
p(s) and z(s) where s is the arclength along the profile (generating curve) of the scatterer. 
The surface of the scatterer is generated by rotating the generating curve in the <p direction 
around the z axis. The surface impedance is allowed to vary as a function of both the 
arclength s and the cylindrical coordinate <p. 

The integral equation formulation leads to vector integro-differential equations 
for the unknown electric surface current density on the scatterer. All three types of 
integral equations, electric field integral equation (EFIE), magnetic field integral equation 
(MFIE) and combined field integral equation (CFIE) are derived. These integral 
equations are then specialized for the case where the scatterer has body of revolution 
symmetry. This specialization transforms the vector integral equations into a pair of 
coupled scalar integral equations in terms of the vector components of the electric surface 
current density. Expressions for the incident fields with arbitrary angle of polarization and 
arbitrary angles of incidence are derived in parallel with the corresponding integral 
equation. Finally, expressions for the far scattered fields and the bistatic radar cross 
sections in terms of the electric surface current density are derived for the geometry under 
consideration. 

In Chapter 3, after a brief discussion of the method of moments to introduce the 
terminology and notations, specific details of geometry and surface impedance 
descriptions are presented. This is followed by a overview of the Hermite expansion 
functions utilized in the method of moments implementation. The discussion covers 
Hermite expansion functions in one-dimension and its extension to two-dimensions. This 
is followed by details of the representation of the surface currents and some practical 
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considerations involving the testing procedure used in this method of moments 
implementation. The discussions presented in this chapter delineate explicitly the 
expressions needed for numerical implementation of the scheme and shed some light on 
the development of the computer code. 

The numerical results obtained from the code developed based on the 
discussions found in Chapters 2 and 3 are presented in Chapter 4. Results validating the 
code have been previously presented in the report for Phase II of this study [Uslenghi et. 
al., 1991b]. The results presented here will concentrate on cases where the surface 
impedance is allowed to have abrupt discontinuities in either the s direction only or in 
both the s and 0 directions (i.e. impedance patch). The data presented was generated on a 
Cray Y-MP with 64 megawords (512 megabytes) of memory. Future plans include the 
parallelization of the code for shared memory parallel processor systems like the Cray Y- 
MP and also for massively parallel processor systems like the Thinking Machines CM5. 



2. ANALYTICAL RESULTS 


2.1 Introduction 

This chapter deals with the derivation of the analytical results for the problem of 
scattering of an electromagnetic plane wave by a scatterer with body of revolution 
symmetry. The material properties of the scatterer are specified through the use of a 
surface impedance boundary condition. General integral equations for a scatterer with 
arbitrary geometry and arbitrary surface impedance are derived. Three integral equations, 
electric field integral equation, magnetic field integral equation and combined field 
integral equation, are presented. These integral equations are then specialized to the 
geometry of interest. Finally, expressions for the far field and bistatic radar cross section 
are presented. 

2.2 Statement of Problem 

Given a scatterer geometry, the material properties of the scatterer and the form 
of the incident electromagnetic fields the problem is to determine the fields which arise 
when the incident fields interact with the scatterer. The geometry of the problem under 
consideration is shown in Figure 2.2.1. The scatterer consists of a body of revolution with 
a surface impedance boundary condition. The surface impedance is allowed to vary 
arbitrarily over the surface of the scatterer. The incident electromagnetic wave is taken to 
be a linearly polarized plane wave with arbitrary angle of polarization and arbitrary 
angles of incidence. Throughout this work, a time harmonic variation of the form e m is 
implicitly assumed. 

The geometry of the scatterer is defined by specifying a generating curve which 
is parameterized in terms of the arclength along the generating curve as shown in 
Figure 2.2.2. To allow for rather complex scatterer geometries the generating curve is 
defined by specifying two profile functions p(r) and z(s) where s is the arclength 
coordinate of the point ( p,z ) on the generating curve and p and z are the usual cylindrical 
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coordinates. The surface of the scatterer is then defined by rotating the generating curve 
about the z axis through an angle of 360°. Since p and z are functions of s, each point on 
the surface of the scatterer is uniquely defined by its s and <f> coordinates. The surface of 
the scatterer can then be easily mapped to a rectangular region in the s-<f> plane. 





Figure 2.2.2 Geometry of generating curve. 
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The direction of propagation of the incident fields is specified by two angles of 
incidence, 6 i and 0, where d i is the angle between the propagation vector, k', and the 
positive z axis and 0, is the angle between the projection of the propagation vector on to 
the x-y plane and the positive x axis. Using 0, and 0,, define the following set of 
orthogonal unit vectors 

x, = pcos^jcosfA^ + ^cosf^jsinfAij-zsinf^) (2.2.1) 

y ; = -psin(A,) + 0cos(A,) (2.2.2) 

z, = psinffl^cosfA^ + ^sin^Osu^A^ + zcos^,) (2.2.3) 

where A, = 0, -0. These unit vectors define a coordinate system in which the incident 
fields propagate along z and have components which lie only in the x ( -y, plane. The 
incident electric field, E',and the incident magnetic field, H', can be written in terms of 
the coordinate system defined by eqs. (2.2.1), (2.2.2) and (2.2.3) as 

E' (r) = [x; cos(v) + y, sin(v)]e”* r (2.2.4) 

H'(r) = — [-x 1 sin(v) + y j cos(v)]e”' k r (2.2.5) 

Zo 

where k =k Q k i =/t Q z., r =pp+zz, Z 0 is the impedance of free space, k 0 is the free 
space wave number and v is the polarization angle. The incident plane wave is said to 
have transverse magnetic polarization (TM polarization) when the incident electric field 
has a component only in the y direction ( v = ±n/2) and transverse electric polarization 
(TE polarization) when the incident magnetic field has a component only in the y_ 
direction (v = 0,/r) 

In general, T] the relative surface impedance can be an arbitrary function of s and 
0. However, some care must be taken so that the definition of r\ does not lead to an 
ambiguous value for the surface impedance at the points s = 0 and s = s ma . For example, 
rj(^,0) = cos(0) leads to an ambiguous surface impedance at s = 0 and s = that is, r\ 
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has an infinite number of values at these two points. However, if T](s,<p) = \i(<t>-n), 
O<0<2tt where u(r) is the unit step function, then 7] has two values at s = 0 and 
s = but this is an acceptable form since it corresponds to a scatterer with a surface 
which is perfectly conducting in the half space y > 0 and has a relative surface impedance 
equal to 1 in the half space y < 0. 


2.3 Derivation of Integral Equations 

In the phasor domain, Maxwell’s equations can be written as follows: 


V x H = J e + jaiD (2.3.1) 

V x E = -J m -jcoB (2.3.2) 

VD=p f (2.3.3) 

V-B = p m (2.3.4) 


where E, H, D, B are the electric field, magnetic field, electric flux density and the 
magnetic flux density respectively, J,, J m , p e , and p m are the electric current density, 
magnetic current density, electric charge density and magnetic charge density 
respectively and a time dependence of e ,m has been assumed. 

A solution of Maxwell’s equations in free space for the electric field, E, and the 
magnetic field, H, in terms of the electric sources, p e and J, and the magnetic sources, 
p m and J m can be written as follows [Elliott, 1981] 

E(r) = -Vd> e (r) - jco\ t (T)-—Vx A m (r) (2.3.5) 

H(r) = -V<D m (r)-;cuA m (r) +— Vx A,(r) (2.3.6) 

/*0 


where 


<I>,(r) 

*.0r) 


J r jG(r,r')p,(r')</V' 

— f G(r.r')p m (r"fdV' 
4^o v 


(2.3.7) 

(2.3.8) 
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A e (r) = -^\ G (r,r')3 t (r')dV' 


A m (r) = ^|G(r,r')J m (r'W' 


(2.3.9) 

(2.3.10) 


£ 0 is the permittivity of free space, p 0 is the permeability of free space, k 0 = a>^]fi 0 £ 0 is 
the free space wave number and G(r,r') is the free space Green’s function which is given 
by 


G(r,r') = G(R) = — ; R = |r - r'| 
R 


(2.3.11) 


The terms <b t , <b m , A t , and A m are referred to as the electric scalar potential, magnetic 
scalar potential, electric vector potential and the magnetic vector potential respectively. 
Additionally, the charge densities p e and p m are related to the current densities J f and 
J m by the following continuity equations 

V • J c (r) + jcop e (r) = 0 (2.3.12) 

V ■ J m (r) + j(op m (r) = 0 (2.3.13) 


The scattering problem under consideration is formulated by invoking the 
equivalence theorem [Harrington, 1961], [Elliott, 1981]. The scatterer is removed and is 
replaced by an equivalent electric surface current density J„, an equivalent electric 
surface charge density p es , an equivalent magnetic surface current density J^, and an 
equivalent magnetic surface charge density p^. The surface currents J„ and and the 
surface charge densities p es and p^ exist on the surface S where S corresponds to the 
surface of the scatterer. The original fields on S and the equivalent surface currents satisfy 
the following boundary conditions 

J =nxH 

es 

J =-nxE 


(2.3.14) 

(2.3.15) 
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where n is the outward pointing surface normal and H and E are the total magnetic field 
and total electric field respectively over the surface S. For the problem under 
consideration the scatterer is taken to be an impenetrable scatterer so the fields are equal 
to zero inside the scatterer. The total electric field and total magnetic field are then given 
by E = E' + E 1 and H = H +H* where E J and H 1 are the scattered electric and magnetic 
fields respectively. 

Since the sources are confined to the surface of the scatterer the potentials 
associated with the equivalent surface current densities and surface charges densities can 
be obtained from eqs. (2.3.7), (2.3.8), (2.3.9) and (2.3.10) by replacing the current 
densities, J f and J m , and the charge densities, p e and p m , by the equivalent surface 
current densities, J Ci and J^, and the equivalent surface charge densities, p es and p ^ 
respectively. Also, the volume integral J dV' is replaced by the surface integral £ dS' . 

Then, the scattered electric and magnetic fields due to the equivalent sources 
can be written as 


E J (r) = 


H'(r) = 


" 7 —VjG(ry)pjT')dS'-^$G(ry)JJr')dS' 

47C£ 0 s 4 71 s 

- — V x | G(r, r') (x')dS' 

471 s 

s ™ s 

+ ^VxjG(r,r')J„(r')iS' 


(2.3.16) 


(2.3.17) 


Interchanging the order of differentiation and integration in eqs. (2.3.16) and (2.3.17), 
making use of the following relations 

Vj G(ry)f(r')dS' = J f(r')VG(r,r')dS ' 

V x J G(r,r')F(r')<tf ' = J VG(r,r') x F(r ')dS' 


(2.3.18) 

(2.3.19) 
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noting that VG = -V'G and using the continuity equations given by eqs. (2.3.12) and 
(2.3.13) to write the charge densities in terms of the surface current densities, the 
scattered electric and magnetic fields can be written as 


E J (r) = ■ JJr')]dS'-^jG(T,r')JJr')dS' 

4 **0 s 5 


f jkp Z 0 


+ -^jVG(r,r')xJ ms (r')dS' 


(2.3.20) 


H r (r) = TT$l V ' G ( r ' r ')P' ■ ^(r')]dS' - £&jG{T,T')J m ,(r')dS' 
4^0 s 4 71 s 

-^jV'G(T,r')x3Jr')dS' 


1/ jko Ip 


(2.3.21) 


where Z 0 = ^fi 0 /£ 0 is the characteristic impedance of free space and T 0 = 1/Z 0 is the 
characteristic admittance of free space. 

If the current densities J w and were known eqs. (2.3.20) and (2.3.21) could 
be used to calculate the scattered electric and magnetic fields. However for a scattering 
problem J ei and J mi are generally not known in advance and so a direct calculation of 
the scattered fields is not possible. Since the incident electric and magnetic fields are 
known it is conceivable that one might try to construct a pair of integral equations for the 
unknown surface current densities by substituting eqs. (2.3.20) and (2.3.21) into the 
electromagnetic boundary conditions given by eqs. (2.3.14) and (2.3.15). However, one 
important piece of information for this problem has not been considered. The material 
properties of the scatterer do not appear in the expression for the either the scattered 
electric field or the scattered magnetic field. They also do not appear in either of the two 
electromagnetic boundary conditions. 

If one were considering the problem of scattering from a dielectric scatterer the 
material properties of the scatterer would appear in the Green’s function used to write the 
potentials for the fields inside the scatterer. For an impedance scatterer the material 
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properties are specified through the use of the surface impedance boundary condition. 
Repeated for convenience, the surface impedance boundary condition is given by 

= E - (E • n)n = 77Z 0 n x H (2.3.22) 

where Z 0 in the impedance of free space, 77 is the relative surface impedance and n is the 
outward pointing surface normal. Taking the cross product of n with eq. (2.3.22) and 
making use of the following relation 

n x E m = n x [e - (E • n)n] = nxE (2.3.23) 

one obtains the following expression 

n x E = T 7 Z 0 n x n x H (2.3.24) 

Making use of the electric and magnetic boundary conditions given by eqs. (2.3.14) and 
(2.3.15) the magnetic surface current density can be written in terms of the electric 
surface current density J„ as follows 

J mj =-7 7 Z 0 nxJ„ (2.3.25) 


Substituting eq. (2.3.25) into eqs. (2.3.20) and (2.3.21) one obtains the following 
expressions for the scattered fields 


E'(r) = 7^|[ v/ C(r,r')][V'. JJr')]dS'-^jG(r,r')JJr')dS’ 

47Tk 0 s 4/r s 

- V ' G( r ’ r ') x fr (r,) "' x 

4 ft s 

H J (r) = --J-j[V'G(r t r')p' ■ (rj(r')n' x J„(r'))]<tf ' 

4 to 0 s 

+ |^G(r,r')[rKr')n' x J„(r ')]dS' 

■^fV'G(r,r')xJ„(r')<tf' 


(2.3.27) 
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The electric field integral equation (EFIE) for the unknown surface current 
density J„ can be derived in the usual manner [Poggio and Miller, 1973] from the 
boundary condition for the total electric field given by eq. (2.3.15). Substituting 
eq. (2.3.26) in eq. (2.3.15), making use of eq. (2.3.25) and taking the limit as the field 
point r approaches r' on the surface S the following electric field integral equation is 
obtained 

77(r)Z 0 n x J„(r) = 2n x E'(r) + x £[V'G(r,r')][V' . JJr')]dS' 

271*0 s 

.Mn X | G(r)r ')j B (r')^ (2.3.28) 

2n J s 

-~nx -f V'G(r,r') x [rj(r')n' x J„(r')]<# ' 

2tt s 

where £ represents the principle value of the surface integral £ . 

In a similar manner, the magnetic field integral equation (MFIE) can be derived 
from the boundary condition for the total magnetic field given by eq. (2.3.14) and the 
expression for the scattered magnetic field given by eq. (2.3.27). The resulting magnetic 
field integral equation is given below. 

3 Jr) = 2h x H'(r) - -£-n x -f [V'G(r,r')|V' • (r7(r')n' x J„(r'))]JS' 

27tk 0 s 

+ x |G(r,r')[7?(r')n' x J„(r'to' (2.3.29) 

2k s 

--^-n x f V'G(r,r') x JJr')dS' 

2.71 s 

Following Mitzner [1968] and Oshiro et. al. [1970], the combined field integral 
equation (CFIE) can be derived from the superposition of the magnetic field boundary 
condition given by eq. (2.3.14) and a modified form of the electric field boundary 
condition given by eq. (2.3.15). This form of the integral equation is also discussed in 
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detail by Jones [1979], Using E,^ = -nxnxE, an alternate form of the electric field 
boundary condition can be written as follows 

nxJ mj = E un (2.3.30) 

Now construct a combined field boundary condition by the superposition of the usual 
magnetic field boundary condition, eq. (2.3.14) and the electric field boundary condition 
given by eq. (2.3.30) as follows 

J„ + S™- (n x ) = n x H + ^E M (2.3.3 1 ) 

A A 

where a CFIE is a scalar constant. The combined field integral equation can then be derived 
by substituting the expressions for the scattered electric and magnetic fields given by 
eqs. (2.3.26) and (2.3.27) into eq. (2.3.31). However, it is more convenient to write the 
combined field integral equation by utilizing the results for the electric and magnetic field 
integral equations. 

Using the notation of Poggio and Miller [1968] the electric field integral 
equation and the magnetic field integral equation can be written in operator notation as 

Ar[J„( r ')] = 2nxE’(r) (2.3.32) 

J w [J.,(r')] = 2nxH i (r) (2.3.33) 

where the electric field integral operator L E and the magnetic field integral operator L„ 
are given by 

MJ„(r')]=n(r)nxJ„(r) 

+ &&L Axicj„(r‘)JS' 

In 3s 

+ ^nx £V'G x [r7(r')n' x J„(r')]dS' 


(2.3.34) 
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and 


£„[J„(r'>] = J„(r) 

+ T^n x iV'C[V' • rj(r')n' x J„(r')]<tt' 

27tk.Q 

ik t (2.3.35) 

-^.nx£r7(r')G[nxJ„(r')]^' 

+^nxf i V'GxJ 0 (r')^' 


The combined field integral equation can then be written in terms of the electric and 
magnetic field integral equations as 


L c [J„(r')] = 2nx 


H'(r) - x E'(r) 


(2.3.36) 


where the combined field integral operator L c is given by 

L c [J„(r')] = L H [J„( r ')]-^nxL f [J„(r')] (2.3.37) 

A 

The electric field, magnetic field and combined field integral equation given by 
eqs. (2.3.32), (2.3.33) and (2.3.36) respectively are valid for scatterers with arbitrary 
shape and arbitrary surface impedance. In the following sections these integral equation 
will be specialized to the case where the shape of the scatterer is a body of revolution. 

2.4 Integral Equations for Body of Revolution 
2.4.1 Preliminary Remarks 

The integral equations obtained in Section 2.3 are valid for an arbitrarily shaped 
scatterer with arbitrary surface impedance. In Sections 2.4.2, 2.4.3 and 2.4.4 these 
integral equations will be simplified for the case of a scatterer with body of revolution 


symmetry. 
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At each point P = (p,<p,z) on the surface of a body of revolution it is possible to 
construct a local right handed coordinate system ( z,n,<f > ) as shown in Figure 2.2.1. The 
unit vectors for the local coordinate system at the point P can be written as 

x = psin(a) + zcos(a) (2.4.1) 

n = pcos(a)-zsin(a) (2.4.2) 

<^ = -xsin(0) + ycos(0) (2.4.3) 

where a is the angle between x and z and p is the usual unit vector from the cylindrical 
coordinate system which is given by 

p = xcos(0) + ysin(0) (2.4.4) 

Also, the unit vectors x and <j) are tangent to the surface at the point P while the unit 
vector n is normal to the surface at the point P. Note, since the generating curve for the 
body of revolution is defined by specifying two profile function p(s ) and z(s ) the angle 
a is also a function of s. 

2.4.2 Electric Field Integral Equation 

In this section the electric field integral equation derived in Section 2.3 for the 
case of an impedance body will be specialized for a scatterer with body of revolution 
symmetry. In the process, the vector integral equation for the unknown vector electric 
surface current density J M will be reduced to a pair of coupled scalar integral equations 
in terms of the vector components of J (1 . 

The electric field integral is specialized by obtaining the components of the 
integral equation which are tangent to the surface of the body of revolution. In order to 
obtain the tangential components it is convenient to rewrite the electric field integral 
equation given by eq.(2.3.28) as follows 


2n x E'(r) = t)(r)Z 0 n x J„(r) - 3f (r) - 3f (r) - 3f (r) 


(2.4.5) 
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where 

( 2 . 4 . 6 ) 

( 2 . 4 . 7 ) 

( 2 . 4 . 8 ) 

and 


3 f ( < r > = -TT ” x f r V'C(r,i r')][‘ V' ■ . ' 

2^fe 0 s 

3 f(r) = -i^n x <f G(r,r')J„(r')<tf ' 

2;r j 

3 f (r) = -|^n x f V'G(r.r') x [tj( r')n' x JJr')]dS' 
2n * 


J„(rO=x 7 r (r') + 07 ,(r') 

= p[sin(aOcos(AV r ( r, ) _s i n (A')<^(> ,/ )] ( 2 . 4 . 9 ) 

+ 4 >[cos(A')./,(r') - sin(a')sin(A' 7 r (r')] + £cos(a' 7 r ( r ) 


where A' = 0'-0. Interchanging the cross product operation with the integration 
operation and making use of the following expression 

V'G = -VG = (r - r ') 1 + y g -*°« = (r-r'jG, (2.4.10) 

R 

where 

G = 1 + jk 0 R e - jko R (2.4.11) 

the expressions for 3f , 3f and 3f can be written as 

3f(r) = (G, (r,r')[n x (r - r')][V' • J„(r ')]dS' (2.4.12) 

2^fc 0 s 

3f(r) = -^&iG(r,r')[nxJJr')]dS' (2.4.13) 

271 s 

3f (r) = -^ f ftfr'K?, (r,r')[n x (r - r') x (n' x J M (r'))]<tf' (2.4. 14) 

Z/T s 

The r and 0 components are obtained by taking the dot product of eq (2.4.5) with x and 
0 respectively. 
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The electric field integral equation given by eq. (2.3.28) can now be written as 
the following pair of coupled scalar integral equations. 

£;(r) = ir 7 (r)Z 0 7 # (r) 

+ TT" ~l Gi(r,r')p'sin(A') 

+ ^ G(r, r') sin(a')sin(A')7 r (r')dS' 

+ | G(r , r') cos(A')7, (r')dS' 

An { 

+ -- -f t?(r')G, (r, r'){[*cos(A') - 2psin 2 (A72)]cos(a') 

47T s 

sin(a')cos(A')}7,(r')dS' 

t7(r , )G 1 (r,r / )^sin(A')7 T (r')^5 / (2.4.15) 


jW . sin(g') 1 

p' ^ p' 




and 


£;(r) = lr7(r)Z 0 7 r (r) 
, jZo 


Ank, 


-J G](r,r'){[zcos(A') - 2psin 2 (A'/2)]sin(a)+ ^cos(a)J 

r , )+ j_aMr2i , 

ds' p' p' dp’ 

- j G(r, r')sin(a)sin(A')7, (r')dS' 

4k s 

ik Z r 

+ 2-4^^G(rT0[sin(a)sin(a')cos(A') + cos(a)cos(a / )]•fr( ^, )^ , 

h71 s 

p(r')Gi(r,r'){p'sin(a)cos(a') -pcos(a)sin(a') 

An J s 

-!;sir\(a)sm(a')}sm(A')J t (r')dS' 

r l( r/ ) G i (r,r'){£ sin(a)cos(A') 


-[z + 2psin 2 (A72)]cos(a)}y r (r , M5' 


(2.4.16) 
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where £ = p' - p, £ = z'-z and dS’ = p'dp'ds' for a body of revolution. 

Finally, from eq. (2.2.4) the rand 0 components of the incident electric field can 
be written as 

E' T (r ) = {[cos(0,)cos(A,)cos(v) -sintA^sintv^sinta) 

- sin(0, ) cos( v) cos(a)}e _/k r 
£J( r) = [cos(0,)sin(A i )cos(v) + cos(A I )sin(v)]e" A ' r 

where 

k' r = • r = £ o [psin(0,)cos(A i ) + zcos(0,)] 

2.4.3 Magnetic Field Integral Equation 

The magnetic field integral equation derived in Section 2.3 for an arbitrarily 
shaped impedance scatterer is now specialized for a scatterer with body of revolution 
symmetry. The method is virtually identical to that used for the electric field integral in 
the the previous section (Section 2.4.2). 

The magnetic field integral equation given by eq. (2.3.29) can be written in a 
more convenient form as follows 

2n x H'(r) = J M (r) - 3"(r)- 3 2 H (r)- 3?(r) (2.4.20) 

where 

3f(r) = --^-n x -f[V'G(r,r')][V' • (rKr')n' x J„(r'))]<tf ' (2.4.21) 

27CK d s 

3?(r) = &-n x J G(r,r')[77(r')n' x J„(r')]<tf ' (2.4.22) 

2tt s 

3?(r) = ~nx -fV'G(r,r') x J„(r')4S ' 


(2.4.17) 

(2.4.18) 

(2.4.19) 


(2.4.23) 
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and J M is given by eq. (2.4.9). Interchanging the cross product operation with the 
integration operation and making use of eq. (2.4.10) the expressions for 3f, 3f and 3f 
can be written as 

3, W(r) = ' fG, (r,r')[n x (r - r')][V' . (rj(r')n' x J„(r'))]<tf ' (2.4.24) 

3 a *(r) = rKr')G(r,rO[nx(«i'x J„(r'))]<tf' (2.4.25) 

3f(r) = ~ -f G, (r, r')[n x (r - r') x J„(r')]<tf' (2.4.26) 

271 s 

The rand <t> components are obtained by taking the dot product of eq (2.4.20) with r and 

A 

<t> respectively. 

The magnetic field integral equation given by eq. (2.3.29) can now be written as 
the following pair of coupled scalar integral equations. 
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tf'(r) = 


4 j - (r) 

- —~r ^G,(r,r'){[^cos(A') - 2psin 2 (A72)]sin(a) + £cos(a)| 

47tk 0 s 

.|^ '>y (a ') + ^|y,(r0d5- 

+ 4^T ^ ^i( r,r ){[^ C 0 S (A / ) “ 2psin 2 (A72)]sin(a)+£cos(a)j 


1 ^ ( l2}y r (r'M5' 


.P' 


j-^7(r')G,(r,r'){[2:cos(A')-2psin 2 (A72)]sin(a) + ^cos(a)} 
4to 0 5 

1 A' p' dp' J 

-^■^( r/ )G( r , r, )[ cos (a) cos (a') + si n (a)si n (a / ) co s(A , )]4( r ')^' 

h7l s 

- ^ <f Ti(r')G(r, r')sin(a) sin( A')J T (r')dS' 

An] 

- - 7 - -f G, (r,r')[p'sin(a)cos(a') - psin(a')cos(a) 

4;r * 


- £sin(a)sin(a')]sin(A')./ T (r , )*£S' 

- ~ j Gj (r,r'){[* cos(a) - £sin(or)]cos(A') 

4/T ^ 

+2p' cos(a')sin 2 (A72)}7 # (r')^' 


(2.4.27) 


and 
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/*;<!•) = |y,(r) 


--^ r fc,(r,r')p'sin(A') 

47Tk 0 s 


,Jn(r-)sin(g-) + a g (r2l ^ 
P' ) 


+ 71 T f C 1 (r,r')p'sin(A')U^^p,( r ')<«' 

4^ J s [p d<f> } 

-—jri(r')G(r,r')sin(a')sin(A')J t (r')dS' 

A7t s 

+ Aj n ( r ')G(r,T')cos(A')J t (T')dS' 
h7C s 

+ j-G, (r, r')4 sin(A')J t (r')dS' 

*\7t s 

+ j- ^G] (r, r'){[^cos(a') - |sin(a')]cos(A') 
47Z s 

-2pcos(a')sin 2 (A'/2)}j r (r')dS' 


(2.4.28) 


where X ~ P' ~P- %-z -z and dS' = p'dp'ds ' for a body of revolution. 

Finally, from eq. (2.2.5) the r and <p components of the incident magnetic field 
can be written as 


//'(r)= -^-{sin(0,)sin(v)cos(a) 

A> 

- sin(cr)[cos(0 j )cos(A, )sin(v) + sin(A i )cos(v)]}e" ;k r 
H‘ t (t)= — [cos(A,)cos(v)-cos(0 1 )sin(A i )sin(v)]e'-' k r 

Z 0 


(2.4.29) 

(2.4.30) 


where k' • r is given by eq. (2.4.19). 

2.4.4 Combined Field Integral Equation 

Since the combined field integral equation for an arbitrarily shaped scattered 
was constructed from the electric and magnetic field integral equations for an arbitrarily 
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shaped scatterer it is reasonable to assume that the combined field integral equation for a 
body of revolution can be derived from the appropriate combination of the electric and 
magnetic field integral equations for a body of revolution. In order to do this we must first 
determine the rand 0 components of the combined field integral equation. 

The r component of the right hand side of eq. (2.3.36) is given by 


2nx H'W-^nxE'fr) =x- 2nx ff(r)-^nxE'(r) 


= 2 tf;(r) + ^E;(r) 


(2.4.31) 


The r component of the left hand side of eq. (2.3.36) is given by 


= J), 


(2.4.32) 


Then, the r components of the combined field integral equation is given by 


+ a OTE £;< r) = -^(i-rU,,)), 


(2.4.33) 


where \(L H [S es ]) r is given by the right hand side of eq. (2.4.28) and — J„ ])^ is 
given by the right hand side of eq. (2.4.16). 

The 0 component of the right hand side of eq. (2.3.36) is given by 


2nx[ H'(r)-— nxE'(r) =0- 2nx H i (r)-^M»xE , (r) 


= 2l-//;(r) + ^£;(r) 


(2.4.34) 
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The <t> component of the left hand side of eq. (2.3.36) is given by 

(2.4.35) 

Then, the <p components of the combined field integral equation is given by 

Z 0 «;(r) - acm E;(r) = (2.4.36) 

where is given by the right hand side of eq. (2.4.27) and is 

given by the right hand side of eq. (2.4.15). 

2.5 Far Field and Radar Cross Section 

2.5.1 Far Scattered Fields 

In the far field the scattered electric and magnetic fields can be written much 
more simply [Elliott, 1981] as 

E'(r) = -jco{\ eT (r)- Z 0 r x A mr (r)) (2.5.1) 

H* (r) = -Jco(A mT (r) + Y 0 r x A er (r )) (2.5.2) 

where A tT and A mT are the transverse (to r) components of the electric and magnetic 
vector potentials and are given by the following expressions 

A er (r) = A,(r) - f(r • A e (r)) = 0A, e (r) + 0A„(r) (2.5.3) 

A mr (r) = A m (r) - f (f • A m (r)) = 04^ (r) + ^ (r) (2.5.4) 

A A 

where f , 0 and <{> are the usual unit vectors form the spherical coordinate system. The 
free space Green’s function can also be simplified in the far field region . When |r| is 
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much larger than |r'| the 1 /R factor, which represents the magnitude of the free space 
Green’s function, can be simply approximated by 1 Jr. For the R factor which appears in 
the complex exponential, the following more exact approximation must be used 


R = |r- r'| «= r - r • r' = r- L 

(2.5.5) 

where, in cylindrical coordinates, 

L - r • r' = p'cos(A')sin(0) + 2'cos(0) 

(2.5.6) 

Thus in the far field the vector potentials can be written as 

A,(r) = -t 2 -— | JJr')e ik ° L dS' 
4 tt r j 

(2.5.5) 

p „-i k o r , 

A.(r) = |a- ij„(r ')e^dS' 

An r \ 

(2.5.6) 


Substituting eqs. (2.5.3) and (2.5.4) into the far field expression for the scattered 
field given by eq. (2.5.1) the 6 and <p components of the scattered electric field can be 
written as 

E’(r) = -ja* ee (r)-jojZ 0 A mt (r) (2.5.7) 

£?(r) = -jCtA^ir) + joiZ 0 A mg (r) (2.5.8) 


Using the far field forms of the vector potentials given by eqs. (2.5.5) and (2.5.6) the 
components of the scattered electric field given by eqs. (2.5.7) and (2.5.8) can be written 
as 


Ei( r) = 


£» = 


h r | 

o~i k o r 


k 0 r 






V 

k 0 r 


-S e (d,<t>) 


■ S,(6,<t> ) 


(2.5.9) 

(2.5.10) 


where 
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^(r') = e-J„(r') 

(2.5.11) 

^(r') = $-J„(r') 

(2.5.12) 

W = 0'J„(r') 

(2.5.13) 

/>') = *• J^(r') 

(2.5.14) 


and S e (8,<l > ) and S t (d,Q) are called the far field scattering (pattern) coefficients. 

Making use of the expression for J„(r') given by eq. (2.4.9) and the expression 
relating to J„ given by eq. (2.3.25) the transverse components of the scattered 
electric field can be written as 

E s e ( r) = - — — ^{[cos(0)sin(a')cos(A')-sin(0)cos(a') 

+ 77(r / )cos(A')] 1 / T (i‘') (2.5.15) 

-[cos(0)sin(A') + 7K r, ) s i n ( a, ) s in(A0]J r # (r , )}e ; *° £ *£S' 

*i 2 

Eg(r ) = ^{[sin(a')sin(A') + T7(r')cos(0)sin(A')]A( r/ ) 

4^" /: 0 r s 

+ [cos(A') + 77(r , )cos(0)sin(a / )cos(A') (2.5.16) 

-r7(^ / )sin(0)cos(a , )]>/ # (^ , )}^’ ; * oi ^ , 


2.5.2 Radar Cross Section 

In the far field, the scattered fields can be written in spherical coordinates as 

E f (r) = 0£» + $E s t (r) = L^-\es e (6,0) + <frS,(0,0)] (2.5. 1 7) 

V 


where S e (6,0 ) and S t (8,<f>) follow from either eqs. (2.5.9) and (2.5.10) or eqs. (2.5.15) 
and (2.5.16). The bistatic radar cross section (differential scattering cross section) is 
defined as [Bowman et. al., 1987] 


<t (6,<{>) = lim4^r 2 

r— 



(2.5.18) 
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Assume a unit amplitude for the incident electric field (i.e. |E'| = 
eq. (2.5.17) the magnitude of the scattered field is given by 


K = E' • E'' = +|S,<9,0)f ] 


Then, the radar cross section can be written as 


o{6,(j>) = lim4 7tr 2 


|S,<ft*)f+|Vft*)f] 4*1, „ a . ,n 


1). From 


(2.5.19) 


(2.5.20) 



3. NUMERICAL METHODOLOGY 


3.1 Introduction 

This section covers the techniques used to find a numerical solution of the 
integral equations derived in Chapter 1. The integral equations are solved using the 
method of moments [Harrington, 1968]. Since the scatterer is a body of revolution the 
surface of the scatterer can be mapped to a rectangular region in the s-0 plane where s is 
the arclength coordinate and 0 is the usual cylindrical coordinate. The rectangular region 
in the s-0 plane is subdivided into many smaller rectangular subregions or patches. On 
each patch the components of the surface current J es (s,0) are expressed in terms of a 
linear combination of two-dimensional Hermite expansion functions. This choice for the 
expansion functions permits the use of delta functions for the testing (weighting) 
functions in the method of moments. 

3.2 Method of Moments 

The electric field integral equation, the magnetic field integral equation and the 
combined field integral equation as by eqs. (2.3.32), (2.3.33) and (2.3.36) all have the 
following general form [Poggio and Miller, 1973] 

L[F(x')] = G(x) (3.2.1) 

where F is an unknown vector function, G is a known vector function and I is a linear 
operator. Construct an approximate solution for F(x) from a linear combination of known 
expansion functions as follows 

N 

F(x) = f(x) = £a„f n (x) (3.2.2) 

where a n are constants to be determined and {f„(x)} is a set of N known expansion 
(basis) functions. 

Define a residual error function 9i(x) as 


27 
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SR(x) = I[f(x')]-G(x) (3.2.3) 

Substituting eq. (3.2.2) for f in eq. (3.2.3) and interchanging the summation and the 
integration the residual error can be written as 

9?(x) = ]£a n L[f„(x')]- G (x) (3.2.4) 

n=l 

Next, define an inner product over S of two vectors x and y as 

(x,y) = JJxy^5 (3.2.5) 

Now choose a set of M weighting (testing) functions (w^) and require that the inner 
product of the weighting functions and the residual error satisfy the following condition. 


( w « = X ( w - ’ (*')]) -< W *.’ G ) = 0 (3.2.6) 

n - 1 

Eq. (3.2.6) can be written in matrix form as 

ZI = V (3.2.7) 

where the elements of Z, I and V are given by 

Z «, = ( w »-Wx')]) (3.2.8) 

h = a n (3.2.9) 

K=(w m , G) (3.2.10) 


The solution to the integral equation can now be found by solving the linear system of 
equations given by eq. (3.2.7) 

The set of expansion functions {!“„} can be completely arbitrary as long as they 
exist in the domain of the operator L. Similarly, the testing functions {wj can be 
completely arbitrary as long as they exist in the range of the operator L. In the problem 
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under consideration evaluating the linear operator applied to an expansion function 
and evaluating the inner product (w m ,L[f,,(x')]) both involve an integration 
over the domain 5 where S is the surface of the scatterer. Therefore, the conditions on the 
expansion and testing functions can be stated more simply as a requirement that they be 
either non-zero over the entire surface of the scatterer (i.e. entire domain S ) or at least 
non-zero over a small portion (i.e subdomain) of the surface of the scatterer. When the 
expansion or testing functions are defined over the entire domain S they are called entire- 
domain basis functions. When the expansion or testing functions are defined over a 
subdomain of S they are called subsectional basis functions. The actual selection of the 
expansion and testing functions is influenced by such factors as the desired numerical 
accuracy, geometry of the problem and expected behavior of the unknown function. 

3.3 Description of Scatterer Geometry and Surface 
Impedance 

In the computer program implementing the solution of eqs. (2.3.46), (2.3.47) 
and (2.3.50) the two profile functions and the relative surface impedance are numerically 
represented by cubic spline interpolation functions [deBoor, 1978]. Since a cubic spline 
interpolation has continuous first and second derivatives it can not be used to describe a 
function which is discontinuous or has discontinuous derivatives. For these cases, it is 
necessary to use more than one spline function to represent the desired function. When a 
function is represented by more than one spline function, the individual spline functions 
are referred to as spline meta-segments. The use of spline meta-segments permits the use 
of a generating curve which produces a surface having geometrical discontinuities (i.e. 
edges) and allows for the use of a discontinuous surface impedance. Note, the spline 
functions used for the generating curve and the spline function used for the surface 
impedance need not be defined over the same set of spline meta-segments. For example, 
for a spherical scatterer with a relative surface impedance rj = 1 for the hemisphere z > 0 
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and r\ = 0 for the hemisphere z < 0 the profile functions p(s ) and z($) would be defined 
using a single meta-segment while the relative surface impedance T] = rj(s) would require 
two meta-segments. 

Finally, the spline description used for the profile functions and the surface 
impedance is particularly useful when it is necessary to evaluate the first and second 
derivatives of the profile functions. For example, it is not necessary to explicitly provide 
the function a(s-) since it only appears in the integral equations in the form of sin(a) and 
cos(a) and these functions are readily available if one makes use of the following: 

cos (a) = — (3.3.1) 

as 

sin(a) = — (3.3.2) 

as 


Similarly, where needed the derivatives of the surface impedance can be easily evaluated. 


3.4 Hermite Expansion Functions 

Suppose that for some unknown function f(x ) we are given the values of /(*), 
dfjdx and possibly higher order derivatives evaluated at a set of points {*„}. The process 
of determining a polynomial interpolation of f(x ) which matches these conditions is 
called Hermitian interpolation [Lancaster and Salkauskas, 1986]. 

Consider the problem of finding an interpolation for f(x ) on a single interval 
[x n ,x H+] ], Given the following values 


/, = /(*„), 


/* + l =/(*„+! X 



(3.4.1) 


the function f(x) can be approximated in the interval [x„,x n+1 ] by the following 


expression 
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/W=/AW+/XW+ Yr+l^/i+1 (■*) fn+\V »+] (x) (3.4.2) 


where <P„(x) and y/„ix ) are cubic Hermite expansion functions given by 

Yr-rY fr-rV 


0»(*) = i 





(3.4.3) 


if x < x < x 


n+1 




V.(*H 


(*-*.) 


(x-x„) 



ifx n . 1 <x<x„ 


i(x„<x< x n+] 


(3.4.4) 


■ J 


and Ax n = x B+1 - x n . This is the simplest polynomial interpolation for f{x) which 
satisfies the boundary conditions given by eq. (3.4.1). 

Using the following properties for <p n (x) and y/„(x) it is easy to show that the 
expansion for f(x) given by eq. (3.4.2) satisfies the boundary conditions specified in 
eq. (3.4.1). 


*(*.)« 1 

^(^-i) = 0«U, + i) = O 
V n (Xn-0 = V'*(X n ) = Vn(X„ +1 ) = 0 

d 

dx 
d¥n 


-1/ r«vvi/ v 

-W =f— ) =0 

,-1 V dx ) n \dx )„+ 1 


(3.4.5) 


" - 


. =1 
dx 

dwA JdVn 

dx J„_, l dx 


= 0 


«+} 


If the values of fix) and df/dx are given for a set of values { x „ } the Hermitian 
interpolation for f(x) is given by 

fix) = £f„(x) = £fn<t>n(x)+ffy„(x) 


(3.4.6) 
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where 


f,w=uw+/xw 


(3.4.7) 


is the Hermite expansion of f(x ) at the point (node) x = x n . Note, the forms of 
eqs. (3.4.2) and (3.4.6) appear to be quite different but the are fundamentally the same. 
Eq. (3.4.2) is simply an expansion for a single interval while eq. (3.4.6) is an expansion 
over multiple intervals. Associating the expansion functions with the nodes makes it 
possible to construct an expansion for a function which is discontinuous or has 
discontinuous derivatives. For example, if f(x) is known to have a discontinuous first 
derivative at x = x n then F n (x) can be written as 


where 


| fA (*) + fn'jVn (*) ^ X < X„ 

1 /AM + fn,,Wr U) ^ * > X m 


(3.4.8) 


fit = 

*-**« ax 


and 


/ n ' r = lim 
*-**« ax 


(3.4.9) 


Similarly, if f(x) is discontinuous at x = x n then F n (x) can be written as 

F ( ) = + iix<x n 

i/».A(^) + />,W if a: > AT#, 

where 

fnj = lim f(x), f n r = lim f{x) 
and f' t and f' r are given by eq. (3.4.9). 


(3.4.10) 


(3.4.11) 


3.5 Hermite Expansion Functions on Rectangles 

For the rectangular region shown in Figure 3.5.1 two-dimensional Hermite 
expansion functions can be obtained from the products of the one-dimensional Hermite 
expansion functions. Define the two-dimensional Hermite expansion functions as follows 



33 


H L( x ^y) = 4>„(x)0 m (x) (3.5.1) 

H L (x,y) = <t>„ (x)y/ m (x) (3.5.2) 

KJx,y) = y/ n (x)(t> m (x) (3.5.3) 

Km(^y) = ¥M)y/Jx) (3.5.4) 


where 0„ and \fr n are the one-dimensional Hermite expansion functions given by 
eqs. (3.4.3) and (3.4.4). The two subscripts n and m indicate that these expansion 
functions are associated with the point (x n ,y m ) (n,m'\h node). The superscripts 1, 2, 3 and 
4 are used to indicate that a particular expansion function is associated with the following 
values 


/ 


3f % 

' dy' ax' dxdy 


evaluated at the n,m'\h node. 



Figure 3.5.1 Domain of two-dimensional Hermite expansion functions. 


Let 


f,j = /(*.y) 


(#} 


UJ 

u dx 


(3.5.5) 




(3.5.6) 
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f j£ 


£ 

dy 


x=x n y=yj 


( & 


of 1 _ of 


dxdy)'. dx dy 


*=*, .y=yj 


(3.5.7) 

(3.5.8) 


Using the expansion functions defined by eqs. (3.5.1), (3.5.2), (3.5.3) and (3.5.4) a 
function f(x,y ) can be approximated in the region x„ <*<*„+,, y m <y<y m+: by the 
expansion 

f(x,y) = F nm (x,y) + F nm+} (x,y) + F n+hm (x,y) + (x,y) 

^ „ , (3.5.9) 

k = 0 (=0 


where 


F u (x,y)= f iJ Hl j (x,y) + 




£ 


Hlj(x,y) 


+(§) H ^ iXty)+ [^k) H - j(x,y) 

\dxJij yfodyjij 


(3.5.10) 


is the Hermite expansion of f(x,y ) at the ij’th node. 

As was the case with the Hermite expansion of a function of one variable, it is 
possible to construct a Hermite expansion for a function of two variable which is 
discontinuous or has discontinuous first derivatives. The only restrictions is that the 
discontinuities must lie along values of constant x or constant y. For example, if df / dx is 
discontinuous along x = x, then the expansions at nodes along x = x i have the following 
form 
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\fij H iA x 'y) + 


§L 

%),j 


Hfj(x,y) 






where 




i.j.l 

f 


% 


fykj 


Hl(x,y) 






.«■# 


L, — dx 

' Jy y=yj 

dxfy) 

1 = lim -f- 


>U.r OX 

V=V. 

[dxdyj 


d 2 f 


JJ 


= lim 

*-*** ax ay 


y=y> 


d 2 f 


= lim 
*-**r ox ay 

>=y> 


if x < x m 


H.j(x,y) ifx>*„ (3.5.11) 


(3.5.12) 

(3.5.13) 


Similarly, if f(x,y ) is discontinuous along x = x, the Hermite expansions at nodes along 
x = Xj have the form 


\<ty) Uj j 




+ -f- 


§L 


2/ \ 


Hl J (x,y) + 


El 

dxd y)i, i . l 


H*j(x,y) 




'§£ 

K^yJi.j.r 


Hlj(x,y) 


+ (% ) Hli(x,y) + 


r EL) 

*%L 


if x < x. 


H*j(x,y) if x> x„ (3.5.14) 


where 


fijj = lim f(x,y) 

J x—*x~ 

y=>j 

fi.j.r = lim f(x,y) 

X-¥X { 

y~ y j 


&kij 

<dy)i,j.r 


lim# 
*->*i ay 

y=>j 

lim — 
*->< dy 

y=y> 


(3.5.15) 

(3.5.16) 
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Similar expressions can be found for the cases when f(x,y) or df / dy are discontinuous 
along y = yj or discontinuous along both x = x t and y = y ; . 

3.6 Expansion of Surface Currents 

Using s-Q coordinates the surface of the scatterer can be mapped to the region 
®- s - 5 niax’ O<0<2n in the s-0 plane. Keep in mind that since we are considering 
bodies of revolution the points (s,0) and (s,2n) in the s-0 plane correspond to the same 
point on the surface of the scatterer. This is a result of the surface of the scatterer being 
periodic in 0 with period 2n. Therefore, we also expect any function defined over the 
surface of the scatterer, such as the components of the surface current, to be periodic in 0 
with period 2 it. Also, note that points along s = 0 represent the same physical point on 
the scatterer. However, the components of the surface current, 7 r and J t are in general 
not constant along 5 = 0 since the unit vectors t and <|> and hence 7 r and 7, are 
functions of 0. These same considerations also apply to points along s = s^. 

The domains(patches) over which the expansion functions for the components 
of the surface are defined are constructed by dividing the region 0 <s< O<0<2n 

as shown in Figure 3.6.1. There are N s - 1 intervals along the j-axis and N t - 1 intervals 
along the 0-axis for a total of (N s - \)(N $ - 1) patches and a total of N S N $ nodes. The 
patches are numbered so that the upper left hand comer of the i./th patch as shown in 
Figure 3.6.1 has the coordinates Cs,,0 ; ). The components of the surface current can be 
written as 

N, N, 

4(^0) = £5 (3.6.1) 

1=1 ;=1 

x, 

1=1 7=1 


(3.6.2) 
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where J Ti j and J $ij are the Hermite expansions of J r and at the point (s,,0 ; ) (i.e. 
i,j ' th node) and have a form similar to either eq. (3.5.10), (3.5.11) or (3.5.14). The 
particular form of the nodal expansion functions will depend upon both the local 
geometry of the scatterer and the continuity of the surface surface at each of the nodes. 


s i s 2 ••• s, ... S s,-1 S S, 



Figure 3.6.1 Expansion function domains for components of the surface 
current. 

Now consider the nodes along the line s = 0 in the s-<j> plane when a(0) = njl . 
For this case the scatterer is smooth and has a well defined normal at the point 5 = 0. The 
point s = 0 on the surface of the scatterer maps to the line s = 0 in the s-</> plane. This is 
due to the fact that the unit vectors X and 0 are not uniquely defined at the point 5 = 0. 
However, since all the points along 5 = 0 represent the same physical point on the 
scatterer it is reasonable to assume that the coefficients for the expansion functions 
associated with nodes along s = 0 are not completely independent. The relationships 
between the coefficients can be determined by selecting one or more coordinate systems 
which have uniquely defined coordinate vectors at s = 0. 
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Construct a new local coordinate system on the surface of the scatterer in the 
neighborhood of s = 0 as shown in Figure 3.6.2. Define a new arclength parameter s' and 
a new tangent vector t* as follows 


• 

II 

cl ^ 

when </> = 0 O 
when 0 = </> o +7t 

(3.6.3) 

r -U 

when = o 
when <f) = 0 o +jt 

(3.6.4) 



Since the generating curve is smooth at s' = 0 the T* component of J tj must be 
continuous at s' = 0. This condition can be stated mathematically as 


.lim T* J „(*»= lim t* J „(s» (3.6.5) 

r -> 0 ' S -40 


Using the definitions of s' and t* eq. (3.6.5) reduces to 


lim JJs,#) = - lim J T (s,0) 

s— j — * 0 * 


(3.6.6) 
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Since a(0) = ;r/2 da/ds * is continuous at s' = 0. Therefore, the derivative of the x 
component of J ej with respect to s at s' = 0 must also be continuous. This condition can 

be written as 


.-‘r I 7 *' ■ = A r ;i 7rj - (J > ) 

(3.6.7) 

For s' < 0 we have x = -x, s' = -s and =-%, and for s' > 0 

/ = s and =%,. Therefore, eq. (3.6.7) reduces to 

we have x = t. 

lim -^-J T (s,<p) = lim ^J T (s,0) 
ds i-+o* ds 

♦=*o+* *-*<, 

(3.6.8) 

The results given in eqs. (3.6.6) and (3.6.8) imply the following: 


d d 

lim — / T (j,0) = -lim — J t (s,<j>) 

00 ►O (J0 

(3.6.9) 

ir ^ W)= J 

(3.6.10) 


For the <f> component of the surface current define a new angular vector 0* as 

follows 

t. f 0 when 0 = 0 O 

<t> = . (3.6.11) 

[ -0 when 0 = 0 O + n 

Like the component of the surface current in the r direction the component in the 0* 
direction must also be continuous at s' = 0. 

Jim 0*-J „(j\0) = lim 0*-J„(i\0) (3.6.12) 

s — >0 s — >0 

♦*»o+* *-*0 

Using the definition for 0 the condition on the continuity of the 0* component can be 


written as 
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lim /,(j,0) = -lim JJs,0) (3.6.13) 

j-kO* i->0‘ T 

#»#o+* *=*o 

Since the geometry of the scatterer is smooth at s = 0 the derivative of the 0’ component 
of J„ with respect to s' should be continuous at s' = 0. 

.lim J Js',0) = lim -£•$*• J Js',0) (3.6.14) 

j as 5 ->o ds 

Since §' = <j) and %• = -% for s* < 0 and $ and %• = %, for s* > 0 eq. (3.6.14) 
simplifies to 

lim 4-y,(r,0)= lim ^-JAs,0) (3.6.15) 

*->CT cfr j-40* cfS 


The results given in eqs. (3.6.13) and (3.6.15) imply that the following is also true. 


d d 

lim — 7 (s,0) = -lim— JAs,0) 

s—*0~ M s-*0* (j(h 


lim 

f=*o + 7r 


dsd0 


JAs,0)= lim 

s — > 0 * 

*=*0 


dsd0 


J.(s,0) 


(3.6.16) 

(3.6.17) 


If aCs,™,) = -n!2 then the results given by eqs. (3.6.6, 8-10, 13, 15-17) also apply in the 
limit s -4 that is, at the point s = s,,^. For nodes along s = 0 (s = s^) whenever 
a( 0) = x/2 (a(s maj ) = -n/2) eqs. (3.6.6), (3.6.8), (3.6.9) and (3.6.10) can be used to 
reduce the number of unknowns per node associated with J T from 4 to 2. Similarly, 
eqs. (3.6.13), (3.6.15), (3.6.16) and (3.6.17) can be used to reduce the number of 
unknowns per node associated with J p from 4 to 2 for all nodes along s = 0 (s = s^) 
whenever a( 0) = n/ 2 (a(s mM ) = -n/2 ). 

If a(0)*?r/2 the scatterer starts with a sharp point. This geometrical 
singularity will results in a possible singularity in J t at s = 0 while J r will remain 
continuous at s = 0. However, since there is not a well defined surface normal at s = 0 
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we can not construct a local coordinate system at s = 0 and hence we can not 
geometrically determine the behavior of either J J T or their derivatives in the limit 
s — * 0 + . These same considerations also apply at s = s ^ when ais^) * -n/2.1n either 
of the two cases eqs. (3.6.6, 8-10, 13, 15-17) can not be used to reduce the number of 
unknowns and 4 unknowns per current component must be used at each node along s = 0 
or * = W 

3.7 Selection of Sampling and Integration Points 

Since Hermite expansion functions belong to the class of functions with C 1 
continuity (i.e. continuous first derivatives) the use of delta functions as the testing 
functions is allowed [Chap. 11, Vichnevetsky, 1981]. From eq. (3.2.13) it can be seen that 
employing delta functions as the testing functions is equivalent to sampling the integral 
equation at some arbitrary point on the surface of the scatterer. Since the domain of the 
testing functions is the surface of the scatterer the sampling points must be located on the 
surface of the scatterer. In addition, since the domains of the expansion functions are a 
subregion of the entire surface, the locations of the sampling points must be distributed 
over the surface of the scatterer in order to to include the effects of all the expansion 
functions. That is, if there are N unknown coefficients then the sampling points must be 
chosen to produce N linearly independent equations in terms of the unknown coefficients. 

Using delta functions as the testing functions, consider the method of moments 
solution to eq. (3.2.7) for the simple case where F and G are the scalar functions f(x ) 
and g(x) respectively. For the interval [x,,**] let the expansion of f(x) have the form 
given by eq. (3.4.6). Assume that at each node x n , n = in the interval [*,,**] the 

function /is continuous so that the nodal expansion functions F n have the form given in 
eq. (3.4.7). Since there are two unknowns per nodal expansion function, the total number 
of unknowns is 2 N. The elements of Z and V given in eqs. (3.2.8) and (3.2.10) are 
determined by sampling L[/(x') ] and g(x) respectively at 2N points in the interval 



42 


Since the expansion function at the n’th node is only defined for the interval 
two sampling points associated with the unknown coefficients from the 
nodal expansion function at x = x n must be located somewhere in the interval 
Since a nodal expansion function contributes to the overall value of f{x ) on both sides of 
the corresponding node, the locations of the sampling points associated with the /fth node 
are chosen so that one sampling point is located in the interval [x„_ lt x n ] while the other is 
located in the interval [*„,*„+,]. The exception to this rule occurs for the first and last 
intervals. Since a solution is desired only for the interval it is not logical to locate 

sampling points outside of this interval. Therefore, the sampling points associated with 
the first and N'th nodes can only be located in the intervals [x,,x 2 ] and 
respectively. The required distribution of sampling points points is shown in Figure 3.7.1. 
If the values of / or df jdx are known at the endpoints or boundary conditions for /and 
df/dx are specified at the endpoints then it may be possible to reduce the number of total 
number of sampling points in the first and last intervals from three to two. 



Figure 3.7.1 Location of sampling points for a continuous function. 

Circles represent nodes; crosses represent sampling points. 


If df/dx is discontinuous at x - x, then a nodal expansion function of the form 
given in eq. (3.4.8) must be used at the Fth node. This introduces an additional unknown 
coefficient for which an additional sampling point must be chosen. This additional 
sampling point must be located at the point x = x, as shown in Figure 3.7.2a. If /( x) is 
discontinuous at x = x, then the nodal expansion function at the Fth node has the form 
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given by eq. (3.4.10). This choice for the nodal expansion function introduces two 
additional unknown coefficients into the problem. While there are four unknown 
coefficients associated with the i’th node two of the coefficients only contribute to /in the 
interval and the other two coefficients only contribute to / in the interval 

1*1’*.+, ]• Therefore, the two of the four sampling points associated with the i’th node are 
located in the interval while the other two sampling points are located in 

[*,’*;+] 1 35 shown in Figure 3.7.2b. 

The distribution of sampling points for the two-dimensional problem is based on 
the distributions found for the one-dimensional problem. For example, consider a two- 
dimensional function f(x,y) where /and df/dy are always continuous in y, df/dx is 
discontinuous in x at x = x t and /is discontinuous in x at x = x r Construct a new one- 
dimensional function /(y) = f(x 0 ,y), where x 0 is a constant. Since /(x,y) and df/dy are 
always continuous in y, / and dfjdy are also continuous in y. Therefore, the distribution 
of sampling point along the y direction should be chosen as shown in Figure 3.7.1. For the 
sampling points in the x direction construct a new function f 2 (x ) = f(x,y 0 ) where y 0 is a 
constant. Since df/dx and f(x,y ) are discontinuous at x = x i and x = x j respectively 
df 2 /dx and f 2 will also be discontinuous at x = x, and x = x ; respectively. Therefore, at 
nodes where /, and df 2 jdx are both continuous choose the sampling points as shown in 
Figure 3.7.1; at nodes where df 2 /dx in discontinuous choose the sampling points as 
shown in Figure 3.7.2a; and at nodes where f 2 is discontinuous choose the sampling 
points as shown in Figure 3.7.2b. 

In principle, the sampling points can be located anywhere in the interval. 
Numerical accuracy is improved if they are located at the points used for Gaussian 
integration [Chap. 11, Vichnevetsky, 1981]. Sampling at the Gaussian points is referred 
to as orthogonal collocation. The drawback in using orthogonal collocation results from 
the complexity it introduces to the choice for the numerical integration rule. 
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(a) (b) 

Figure 3.7.2 Location of sampling points for discontinuous functions, (a) 
df/dx is discontinuous at x = ^/b) f(x ) is discontinuous 
at * = *;. Circles represent nodes; crosses represent 
sampling points. 

The free space Green's function, G, which appears in the integral equations is 
singular whenever the observation point (i.e. sampling point) and integration point 
coincide. The free space Green's function appears in three forms in the integral equation, 
as GF the scalar product of G and a vector function F, as /VG the product of the scalar 
function/ and VG and as VG x F the vector cross product of VG and a vector function 
F. Analytically, if F is non-singular then the integral of GF over a surface is always 
convergent [Kellog, 1953]. However, there are convergence problems associated the 
numerical integration of a function such as GF even when it has been established 
analytically to be integrate. Analytically, the integrals of /VG and VGxF over a 
surface are, in general, not convergent. However, the problems associated with these two 
singular integrals have been accounted for through the use of the concept of principle 
value integrals. 

With this in mind, there are two general possibilities for the choice of the 
locations of the sampling points and integration points. The sampling points and 
integration points can be chosen so that there is a possibility that they will be coincident 
or they can be chosen so that they are never coincident. If the sampling and integration 
points are allowed to be coincident then some provisions must be made in the numerical 
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integration routine to detect the case where the two points are coincident. For the integrals 
of GF the integration in the vicinity of the singularity is either replaced by a numerical 
integration rule specifically designed to handle this type of singularity or evaluated using 
an approximate analytical expression for the integral. For the integrals of /VG and 
VG x F , since the effects of their singularities have already been accounted for through 
the use of the principle value integrals, the integrands can simply be set to zero whenever 
the integration and sampling points are coincident. A simpler procedure is to select the 
sampling and integration points so that they are never coincident. Then, the numerical 
integration does not have to check for the possibility of a singular integrand. There is 
some loss of numerical accuracy if the singularity of GF is not considered explicitly. 
However, this effect can be minimized by selecting the integration points so that they are 
symmetrically placed around the sampling points. 

It is relatively easy to construct a numerical integration rule which will avoid 
selecting integration points which are coincident with the sampling points. Integration 
over a surface patch is performed by dividing the patch into a number of square 
subpatches. The center point of each patch is referred to as an integration point. On each 
integration subpatch the value of the integral is approximated as the product of the area of 
the subpatch times the value of the integrand evaluated at the integration point for the 
subpatch. Let N mx be the number of integration intervals along the x and y axes. 
Assuming that the integration is being performed over the normalized patch -1 < x < 1 
and -1 < y < 1 the numerical integration rule just described corresponds to selecting the 
two sets of weighting coefficients {w^} and {w y/ } for a two-dimensional numerical 
integration rule as 


where k = 1 N m , and selecting the two sets of integration points {jtJ and {>,} as 
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2k -\ , 

x k =y k = 1 

AL 


(3.7.2) 


where /: = . If the value of N m is restricted to a multiple of 4 and the possible x 

and y coordinates of the sampling points in the normalized patch are restricted to the 
values -1/2, 0 and 1/2 then the integration points and the sampling points will never be 
coincident. Note, using this scheme for the selection of the sampling and integration 
points does not permit the construction of a solution which is discontinuous at both ends 
of an interval. 



4. NUMERICAL RESULTS 


4.1 Organization of Computer Code 

The computer code is divided into five distinct programs. Depending on the 
choice of the geometry, a total of four programs will actually be used. The programs are: 

1. CNSPH3 

2. CNCL3 

3. PREMOM3 

4. MOM3 

5. PSTMOM3 

The first two programs, CNSPH3 and CNCL3, are used to define the geometry for the 
cone-sphere and cone-cylinder scatterers respectively. The cross sections for both of these 
geometries are shown in Figures 4.1.1 and 4.1.2. Note that a spherical scatterer is a 
special case of both these geometries. Also, the cone-cylinder-sphere geometry can be 
generated as a particular case of the cone-cylinder geometry. The output from each of the 
two geometry programs is the cubic spline description for the generating curve and a 
cubic spline description for the surface impedance. 

The program PREMOM3 reads the output file generated by one of the geometry 
programs, generates a surface patch model of the scatterer and then generates any 
remaining information needed to solve any of the three integral equations. PREMOM3 
also generates some additional information used in the calculation of the far field 
scattering coefficients and the bistatic radar cross sections. The function of PREMOM3 is 
very similar to the mesh generating function found in commercial finite element 
programs. However, as a mesh generator PREMOM3 is limited to generating a single 
type of curved surface element only for geometries with body of revolution symmetry. 
PREMOM3 generates two output files. The first file is an unformatted file containing the 
surface patch model used by both MOM3 and PSTMOM3. This file also contains 
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information used only by PSTMOM3 in the calculation of the radar cross sections. The 
second file is a formatted file containing a Fortran PARAMETER statement which 
defines the absolute minimum array sizes needed by MOM3 and PSTMOM3. Using the 
PARAMETER statement generated for a particular problem in MOM3 and PSTMOM3 
will always result in the use of the absolute minimum amount of memory. 

As it’s name suggests, MOM3 implements a method of moments solution of any 
one of the three integral equations presented in Chapter 2 using Hermite expansion 
functions with delta function testing. MOM3 reads the surface patch model generated by 
PREMOM3, generates the impedance matrix Z for the desired integral equation and then 
solves the matrix equation given by eq. (3.2.7) for both TE and TM polarizations. Finally, 
MOM3 generates an unformatted file containing the coefficient vectors for the surface 
current belonging to the two orthogonal polarizations and a copy of the surface patch 
model generated by PREMOM3. 

PSTMOM3 reads the file generated by MOM3 and computes, for both 
polarizations, the surface current at all integration points used in the evaluation of the 
scattering coefficients. The far field scattering coefficients are then calculated for a 
specified number of values for 6 and 0. The scattering coefficients are then used to 
calculate the bistatic radar cross section at the same values of 8 and 0. Finally, the results 
of all these calculations, the surface currents, the scattering coefficients and the radar 
cross sections are written to a formatted file. 

An attempt has been made to write the programs so that they adhere as closely 
as possible to the Fortran-77 standard. The programs have been successfully compiled 
without modification by a large number of Fortran compilers ranging from Watfor-77 and 
Microsoft Fortran under MS-DOS to the cf77 compiler system under UNICOS. 
Additionally, some time has been spent improving the computational efficiency of the 
routine in MOM3 which generates the elements of Z (i.e. matrix fill routine). This effort 
has been aimed mainly at writing a matrix fill routine which can be easily ported to vector 
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supercomputers like the Cray Y-MP. Future plans include the parallelization of the code 
for shared memory parallel processor systems like the Cray Y-MP and also for massively 
parallel processor systems like the Thinking Machines CM5. 



Figure 4.1.1 Cross section of cone-sphere geometry 



Figure 4. 1 .2 Cross section of cone-cylinder geometry 
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4.2 Sample Results 

In this section the bistatic radar cross sections for two different cone-cylinder 
scatterers with surface impedances which vary along both the s and </> directions will be 
presented. All results have been obtained by solving the magnetic field integral equation. 
Additionally, all solutions have used 4 integration points per interval in both the s and 0 
directions for a total of 16 integration points per surface patch. 

The first three examples are for a cone-cylinder (see Figure 4.1.2) with 
normalized dimensions k 0 h cy] =3.0, k 0 r cyl =0.5, =21°, k 0 r tip =0.1875, k 0 r n =2.0 

and k 0 r J2 =0.0. This geometry will be referred to as cone-cylinder 1 . The normalized 
arclength for this scattered s generating curve is k^^ = 4.6242 . Results for three 
different variations in the surface impedance of cone-cylinder 1 will be presented. For all 
three examples the surface patches are generated using N s = 13 nodes in the s direction 
and N t =9 nodes in the <j> direction. For the first example, the relative surface impedance 
has a constant value of 77 = 0 (i.e. perfectly conducting). The bistatic radar cross sections 
for this constant surface impedance are shown in Figures 4.2.3 and 4.2.4. The radar cross 
sections are for the case of oblique incidence with angles of incidence 0 , = 0 ° and 
6 t = 45° and both TE and TM polarizations. 

For the second and third examples using cone-cylinder 1 the relative surface 
impedance in piecewise constant with two step discontinuities, one at k 0 s = 0.7619 and 
another at k 0 s = 1.4949 . The surface impedance for the second example is given by 


fo 


r?(V) = 1 


1 

0 


0.0 <V <0.7619 
0.7619 < k 0 s < 1.4949 
1.4949 < k 0 s < 4.6242 


(4.2.1) 


while the surface impedance for the third example is given by 
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n 


ri(k 0 s) = \ 


2 

1 


0.0 < Jt 0 j <0.7619 
0.7619 < k 0 s < 1.4949 
1.4949 < k 0 s < 4.6242 


(4.2.2) 


A graphical picture of the surface impedances given by eqs. (4.2.1) and (4.2.2) are shown 
in Figures 4.2.1 and 4.2.2 respectively. It can be seen from Figures 4.2.1 and 4.2.2 that 
points with a normalized arclength coordinate k 0 s in the range 0.7619 < k 0 s < 1.4949 are 
located on the rounded joint between the cone and the cylinder. The bistatic radar cross 
sections for cone-cylinder 1 with a relative surface impedance given by eq. (4.2. 1 ) are 
shown in Figures 4.2.5 and 4.2.6. The bistatic radar cross sections for cone-cylinder 1 
with a relative surface impedance given by eq. (4.2.2) are shown in Figures 4.2.7 and 
4.2.8. The results presented are for TE and TM polarizations and angles of incidence 
0,=O° and 0,=45°. 

The second cone-cylinder geometry has normalized dimensions k 0 h cyi = 4, 
k 0 r cy] = 1 , a cone =26.565°, k 0 r Up = k a r ]2 =0.25 and k 0 r n =0.75. This geometry will be 
referred to as cone-cylinder 2. The normalized height for cone-cylinder 2 is h = 5.6910 
and the normalized arclength of its generating curve is s ^ =6.8992. Two different 
examples for the surface impedance will be used with this geometry. The first example 
uses a constant relative surface impedance r\ = 1 (i.e. perfectly absorbing). For this case 
the surface patches have been generated using N s = 15 and N* = 13 which results in a 
total of 168 surface patches and a total of 1584 unknowns in the method of moments 
program. Results for the bistatic radar cross sections for this scatterer due to both TE and 
TM polarized incident fields with angles of incidence <t> t =270° and 0, =0°, 45°, 90°, 
135° and 180° are shown in Figures 4.2.10-19. The reason for choosing <f> i = 270° is so 
that the results for cone-cylinder 2 with T] = 1 can be directly compared with the results 
obtained for the next example. 

The second example using cone-cylinder 2 uses a relative surface impedance 
77 = 1 for the entire surface except for a rectangular region where the relative surface 
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impedance is r\ — 0. The surface impedance for this example is given by the following 
expression 


n(V.0)=i-f,(* o *)W 


where the functions f, and f # are given by 


fo 


f x (V) = i 


1 

0 


0 < IcqS < 3.1835 
3.1835 < k 0 s < 4.6835 
4.6835 < k 0 s < 6.8992 


f,(0) = 


0 < 0 < n/4 
nj4 < 0 < 3?r/4 
2>n/4 < 0 < 2 n 


(4.2.3) 


(4.2.4) 


(4.2.5) 


A side view of cone-cylinder 2 with the surface impedance given by eq. (4.2.3) is shown 
in Figure 4.2.9. Results for the bistatic radar cross sections for cone-cylinder 2 with the 
relative surface impedance given by eq. (4.2.3) due to both TE and TM polarized incident 
fields with angles of incidence 0, = 270° and 0, = 0°, 45°, 90°, 135° and 180° are shown 
in Figures (4.2.20-29). The reason for selecting 0_ = 270° is so that an observer looking at 
the scatterer along the propagation vector k‘, where k‘ = k 0 k i = fc 0 z and z, is given by 
eq. (2.2.3), will see the impedance patch. 



53 


k 0 s= 1.4949 



k 0 s 


4.6242 


Figure 4.2.1 Graphical representation of relative surface impedance 
given by eq. (4.2.1) 


k n s = 1.4949 



k,s = 4.6242 


Figure 4.2.2 Graphical representation of relative surface impedance 
given by eq. (4.2.2) 



RCS (dB) 
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Figure 4.2.3 Bistatic radar cross section for perfectly conducting cone-cylinder: 
TE polarization, Q t = 0° and 0,=45° - * 0 /i , =3.0, /t 0 ri=0.5, 

« cone = 21°, * 0 r, up = 0.1875, k 0 r n = 2.0, k 0 r J2 = 0.0, N t = 13, and N $ = 9. 
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<r> 



Figure 4.2.4 Bistatic radar cross section for perfectly conducting cone-cylinder: 
TM polarization, 0,=O° and 0 ; =45° - k 0 h , = 3.0, k 0 r.= 0.5, 

a OT£ = 21°, Vu P =0.1875, k 0 r n =2.0, k 0 r J2 =0.0, N t = 13, and N $ = 9,. 


RCS (dB) 
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Figure 4.2.5 Bistatic radar cross section for cone-cylinder with piecewise constant 
surface impedance: TE polarization, 0 ( =0° and d i =45° - k 0 h , = 3.0, 
Vcyi=0-5, a cme =21°, Vup = 0.1875, k Q r n = 2.0, k Q r J2 =0.0, N, = 13 
and N t = 9 with relative surface impedance given by eq. (4.2.1). 


RCS (dB) 
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Figure 4.2.6 Bistatic radar cross section for cone-cylinder with piecewise constant 
surface impedance: TM polarization, 0, = 0° and 6 \ = 45° - Lh . = 3.0, 
^o^cyi — 0*5, cc cone — 21 , ko r w = 0.1875, k 0 rj i = 2.0, £ 0 r, 2 = 0.0, N t = 13, 
and N t = 9 with relative surface impedance given by eq. (4.2.1). 


RCS (dB) 
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Figure 4.2.7 Bistatic radar cross section for cone-cylinder with piecewise constant 
surface impedance: TE polarization, 0, =0° and 6 =45° - k 0 h , = 3.0, 
k 0 r cy] = 0.5, a cone = 21°, k 0 r = 0.1875, k 0 r n =2.0, k 0 r J2 = 0.0, ft, = 13, 
and N'=9 with relative surface impedance given by eq. (4.2.2). 


RCS (dB) 
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Figure 4.2.8 Bistatic radar cross section for cone-cylinder with piecewise constant 
surface impedance: TM polarization, <p t =0° and 6, = 45° - k 0 h , = 3.0, 
Vcyi = 0.5, a cone =21°, k 0 r = 0.1875, k Q r n = 2.0, k 0 r J2 = 0.0, ^, = 13, 
N t =9, <f>' = 0°, and 6, = 45° with relative surface impedance given by 
eq. (4.2.2). 


3.1835 ^ „ s =4.6835 
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RCS (dB) 
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Figure 4.2. 10 Bistatic radar cross section for perfectly absorbing cone-cylinder: 
TE polarization, 0, =270° and 0,=O° - k 0 h ,=4, k 0 r cvl =1, 

a*™ = 26.565°, V„ P = k 0 r n = 0.25, Vji = 0-25, N t = 15 and TV, = 13. 


RCS (dB) 
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Figure 4.2. 1 1 Bistatic radar cross section for perfectly absorbing cone-cylinder: 
TM polarization, 0,=27O° and 0,=O° - *<A y i = 4> *o r cyi = l> 

a*** = 26.565°, k 0 r^ = k 0 r }2 = 0.25, k 0 r n = 0.75, N t = 15 and N f = 13. 



RCS (dB) 
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Figure 4.2. 12 Bistatic radar cross section for perfectly absorbing cone-cylinder: 
TE polarization, 0,=27O° and 0,=45° - k 0 h ,=4, k 0 r cy] =1, 

«core = 26.565°, * 0 r tip = k 0 r J2 = 0.25, Vj, = 0.75, N t = 15 and N f = 13. 


RCS (dB) 
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Figure 4.2. 13 Bistatic radar cross section for perfectly absorbing cone-cylinder: 
TM polarization, 0,=27O° and 0,=45° - k 0 h } = 4, k 0 r. =1, 

0C ccnc = 26.565°, £ 0 r Up - k 0 r 12 = 0.25, k 0 r n = 0.75, N t = 15 and N, = 13. 



RCS (dB) 
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Figure 4.2. 14 Bistatic radar cross section for perfectly absorbing cone-cylinder: 
TE polarization, 0. =270° and 0,=9O° - * 0 /t , = 4, k 0 r cyl = 1, 

« conI = 26.565°, Vu P = /: 0 r J2 = 0.25, k 0 r n = 0.75, JV, = 15 and N, = 13. 
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Figure 4.2. 15 Bistatic radar cross section for perfectly absorbing cone-cylinder: 
TM polarization, 0, =27O° and 0, =90° - k. 0 h x =4, k 0 r. = 1 , 

a cont = 26.565°, Vu P = *o r J 2 = 0.25, Vji ~ 0.75. N, = 15 and N t = 13. 


RC5 (dB) 
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Figure 4.2. 16 Bistatic radar cross section for perfectly absorbing cone-cylinder: 
TE polarization, 0,=27O° and 0,=135° - k 0 h ,=4, k 0 r , ■ , = 1, 
& cane = 26.565°, k 0 r Up = * 0 r J2 = 0.25, Vu = 0.75, N t =15 and N, = 13. 


RC5 (dB) 
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Figure 4.2. 17 Bistatic radar cross section for perfectly absorbing cone-cylinder: 
TM polarization, 0,=27O° and 0,=135° - k 0 h ,=4, l^r. =1, 
Ofecne = 26.565°, k 0 r Up = k 0 r l2 = 0.25, k Q r n = 0.75, N s = 15 and N $ = 13. 


RCS (dB) 
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Figure 4.2. 18 Bistatic radar cross section for perfectly absorbing cone-cylinder: 
TE polarization, 0 1 =27O° and ^=180° - ^ 0 /i, = 4, b,=l, 

a*™ = 26.565°, Vup = Vi 2 = 0.25, Vj. = 0.75, = 15 and N, =13. 
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Figure 4.2. 19 Bistatic radar cross section for perfectly absorbing cone-cylinder: 
TM polarization, 0 1 =27O° and #, =180° - k 0 h , = 4, k 0 r. = 1 , 
= 26.565°, k 0 r^ p = k 0 r n = 0.25, k 0 r n =0.75, N s = 15 and N $ = 13. 


RCS (dB) 
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Figure 4.2.20 Bistatic radar cross section for cone-cylinder with impedance patch: 
TE polarization, 0,=27O° and 0,=O° - k 0 h^ = 4, * 0 r, = 1, 

«cone = 26.565°, k 0 r u = k 0 r S2 = 0.25, k 0 r }] = 0.75, N t = 16 and N, =11; 
with r\ given by eq. (4.2.3). 
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Figure 4.2.21 Bistatic radar cross section for cone-cylinder with impedance patch: 
TM polarization, 0,=27O° and 6,= 0° - k 0 h cy] = 4, k 0 r, =1, 
ttcon. = 26.565°, k 0 r x , = k 0 r J2 = 0.25, k Q r sl = 0.75, N, = 16 and N f = 11; 
with 71 given by eq. (4.2.3). 
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Figure 4.2.22 Bistatic radar cross section for cone-cylinder with impedance patch: 
TE polarization, <p t = 270° and 0,=45° - k 0 h ,=4, * 0 r cyI =l, 

"core = 26.565°, Vuo = Vj 2 = 0.25, Vji= 0.75, W, = 16 and W, = 11; 
with r\ given by eq. (4.2.3). 


RCS (dB) 
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Figure 4.2.23 Bistatic radar cross section for cone-cylinder with impedance patch: 
TM polarization, 0,=27O° and 0,=45° - k 0 h ^=4, jfc 0 r , = 1, 
a «ne = 26.565°, k 0 r. = k 0 r 12 = 0.25, k 0 r , , = 0.75, N s = 16 and N f = 11; 
with r\ given by eq. (4.2.3). 
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Figure 4.2.24 Bistatic radar cross section for cone-cylinder with impedance patch: 
TE polarization, 0, =270° and 0, = 90° - k 0 h , =4, ^, = 1, 

«con e = 26.565°, k 0 r^ =k 0 r n =0.25, k 0 r n =0.75, N t = 16 and N # =ll; 
with T] given by eq. (4.2.3). 


RCS (dB) 
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Figure 4.2.25 Bistatic radar cross section for cone-cylinder with impedance patch: 
TM polarization, 0,=27O° and 0,=9O° - k 0 h , = 4, k 0 r x = 1, 

a«n. = 26.565°, Vu P = *062 = °- 25 > Vji= 0 - 75 . N s - 16 and N, = 11; 
with 77 given by eq. (4.2.3). 
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Figure 4.2.26 Bistatic radar cross section for cone-cylinder with impedance patch: 
TE polarization, 0, = 27O° and #, =135° - i 0 L=4, Ay , :=1, 

= 26.565°, k Q r u = k Q r I2 = 0.25, k 0 r n =0.75, N s = 16 and N, = 11; 
with Tj given by eq. (4.2.3). 


RCS (dB) 
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Figure 4.2.27 Bistatic radar cross section for cone-cylinder with impedance patch: 
TM polarization, = 270° and 0^=135° - k Q h x = 4, £ 0 r.=l, 

a*™ = 26.565°, k 0 r u = k 0 r n = 0.25, ^=0.75, N, = 16 and N, = 11; 
with T) given by eq. (4.2.3). 


RCS (dB) 
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Figure 4.2.28 Bistatic radar cross section for cone-cylinder with impedance patch: 
TE polarization, 0,=27O° and 0, =180° - * 0 /r 1 =4, Jfc 0 r cyl = 1, 

a c = 26.565°, k 0 r u = k 0 r n = 0.25, k 0 r n = 0.75, N s = 16 and N, = 11; 
with r\ given by eq. (4.2.3). 


RCS (dB) 
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Figure 4.2.29 Bistatic radar cross section for cone-cylinder with impedance patch: 
TM polarization, 0,=27O° and 0,=18O° - k 0 h t ,=4, k 0 r. = 1, 
a cone = 26.565°, k Q r u = k Q r i2 = 0.25, k 0 r n = 0.75, N t = 16 and N, = 1 1; 
with T] given by eq. (4.2.3). 
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